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The thermodynamical properties of tlie Reissner-Nordstrom-anti-de Sitter black hole in the grand 
canonical ensemble are investigated using York's formalism. The black hole is enclosed in a cavity 
with finite radius where the temperature and electrostatic potential are fixed. The boundary condi- 
, tions allow one to compute the relevant thermodynamical quantities, e.g. thermal energy, entropy 

• and charge. The stability conditions imply that there are thermodynamically stable black hole so- 

I lutions, under certain conditions. By taking the boundary to infinity, and leaving the event horizon 

04 . and charge of the black hole fixed, one rederives the Hawking-Page action and Hawking-Page specific 

5^ ' heat. Instantons with negative heat capacity are also found. 
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I. INTRODUCTION 

m 

' The path- integral approach to the thermodynamics of black holes was originally developed by Hawking et al. . 
, In this approach the thermodynamical partition function is computed from the path-integral in the saddle-point 
approximation, thus obtaining the thermodynamical laws for black holes. 

In the path-integral approach we can use the three different ensembles: microcanonical, canonical and grand 
canonical. Due to difficulties related to stability of the black hole in the canonical ensemble, the microcanonical 
\ ensemble was originally considered However, further developments by York et al. |^-^ allowed to define the 

canonical ensemble. Effectively, by carefully defining the boundary conditions, one can obtain the partition function 
of a black hole in thermodynamical equilibrium. This approach was further developed to include other ensembles [0, 
and to study charged black holes in the grand canonical ensemble [|ll| and black holes in asymptotically anti-de Sitter 
spacetimes [p^p|Jl^. This approach was also applied to black holes in two [Q and three |g2| dimensions. 

In York's formalism the black hole is enclosed in a cavity with a finite radius. The boundary conditions are 
defined according to the thermodynamical ensemble under study. Given the boundary conditions and imposing the 
appropriate constraints, one can compute a reduced action suitable for doing black hole thermodynamics pl] , p^ . 
Evaluating this reduced action at its stable stationary point one obtains the corresponding classical action, which is 
' related to a thermodynamical potential. In the canonical ensemble this thermodynamical potential corresponds to the 
Helmholtz free energy, while for the grand canonical ensemble the thermodynamical potential is the grand canonical 
potential From the thermodynamical potential one can compute all the relevant thermodynamical quantities 

and relations [16[. 

Some controversy has appeared related to the boundary conditions chosen in this formalism ||l^,|l^,^. More 



precisely, Hawking and Page [y_7 18| fix the Hawking temperature of the black hole (i.e. the ternpcrature defined so 
that the respective Euclidean metric has no conical singularity at the horizon), while York et al. ]^,|ri|,Q fix the local 
temperature at a finite radius, where the boundary conditions are defined. For asymptotically fiat spacetimes the two 
formalism coincide, because at infinity the local temperature is equal to the Hawking temperature. On the contrary, 
for asymptotically anti-de Sitter spacetimes the two procedures disagree, since the local temperature is redshifted to 
zero at infinity and is equal to the Hawking temperature only in the region where spacetime has a flat metric. Louko 
and Winters-Hilt jl^] have studied the thermodynamics of the Reissner-Nordstrom-anti-dc Sitter black hole fixing a 
renormalized temperature at infinity that corresponds to the same procedure used in [ pT|JT^ ]. In this paper we have 
chosen to follow York's formalism y,|ll],^ and study the thermodynamics of the Reissner-Nordstrom-anti-de Sitter 
black hole fixing the local temperature at finite radius. 

We find that the two procedures give some identical results, e.g., in both procedures the Hawking-Bekenstein 
formula for the entropy ]l9| , p0t is obtained. In addition, by that taking the boundary to infinity, and leaving the event 
horizon and charge of the black hole fixed, we rederive the Hawking-Page action and Hawking-Page specific heat from 
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York's formalism. However, the value for the energy at infinity differs depending on which procedures one uses. In 
1^ it was found that the energy at infinity is equal to the mass of the black hole, a result that does not hold here. 
These results conform with the similarities and differences found for the Schwarzschild-anti-de Sitter black hole in 

In section O we briefly introduce York's formalism. In section III we compute the reduced action for the Reissner- 
Nordstrom-anti-de Sitter black hole and evaluate its thermodynamical quantities. In section IV we analyze the black 
hole solutions. In section ^ we study the local and global stability of these solutions. The limit where the boundary 
is taken to infinity mentioned above is studied in section VI . Finally some special cases are briefly referred in section 



II. THE ACTION 

The Euclidean Einstein-Maxwell Q is given by 

7 = --^ / d^x^{R-2A) + ^ f dWhK - f d^x^F^.F^''' - Js^btr , (1) 

167r Jj^ Sir Jqm 167r 

where M. is a, compact region with boundary dAA, R is the scalar curvature, A the cosmological constant, g the 
determinant of the Euclidean metrics, K the trace of the extrinsic curvature of the boundary dM, h is the determinant 
of the Euclidean induced metrics on the boundary, F^i, — d^A^, — d^Afj_ is the Faraday tensor and /subtr is an arbitrary 
term that can be used to define the zero of the energy as will be seen later. 

In order to set the nomenclature we follow pT| ] in this section. We consider a spherical symmetric static metric of 
the form ^ 

ds^ = b^dr^ + a^dy^ + r^dO^ , (2) 

where a, h and r are functions only of the radial coordinate y g [0, 1]. The Euclidean time t has period 27r. The event 
horizon, given by y = 0, has radius r+ = r(0) and area = Anr'^. The boundary is given by y = 1 and at this 
boundary the thermodynamical variables defining the ensemble are fixed. The boundary is a two-sphere with area 
Ab = 47rr^, where rs = We will consider the grand canonical ensemble, where heat and charge can flow in and 
out through the boundary to maintain a constant temperature T = T{rB) and electrostatic potential (j) = 0(^-b) at 
the boundary. 



We impose a black hole topology to the metric (0), by using the conditions, 5(0) = 0, ^ =1 and f ^ ) 



see 0- ^ ^ 

Evaluating the action (|l|) for the metric (^, 




ab + Kabr^ \ - - dr 



-0, 

y=o 
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" y=o 

7subtr • (3) 

In order to obtain the reduced action one uses the proper constraints. For the gravitational part of the action Ig 
given in (||), the constraint used is the Hamiltonian constraint lll|Jl5[] 

C'r + Ag^r = SttT^^ . (4) 

In addition, for the matter fields part of the action we use Maxwell equations F^'';^ = 0. 

The thermodynamical quantities and relations are obtained from the "classical action" / (defined as the reduced 
action evaluated at its locally stable stationary points) using the well known relation between the "classical action" 
and the thermodynamical potential 

i = PF . (5) 

Here F is the grand canonical potential since we are considering the grand canonical ensemble. All the thermody- 
namical quantities can be obtained from F using the classical thermodynamical relations (see for example [|l6| ). 
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III. THE reissner-nordstr6m-anti-de sitter black hole 



The Reissner-Nordstrom-anti-de Sitter black hole in the grand canonical ensemble is obtained using a negative 
cosmological constant A and the boundary conditions T = T{Rb) and (j) = 4'{''^b), where is the boundary radius 
of the spherical cavity, T the temperature at the boundary and 4> the electrostatic potential difference between the 
horizon and the boundary. Instead of T we can also use its inverse (3. 

The reduced action for Reissner-Nordstrom-anti-de Sitter black hole is obtained from the Euclidean Einstein- 
Hilbert-Maxwell action / given in (^. For simplicity we split the action in two terms I = Ig + Im, where Ig is the 
gravitational term and the matter field term. To obtain the reduced action we use the Hamiltonian constraint (^) 
and the Maxwell equations. 

The evaluation of Im is identical to the case A = and can therefore be found in lll^l, 



■J3e(j), 



(6) 



where e is the electrical charge of the black hole and (p is the difference of potential between y = and y — 1. To 
evaluate the gravitational term Ig 



we use as mentioned above the constraint 
G\+Ag\ = 8ttT\ . 
The component of the Einstein tensor G^r for the metric (|^) is 



_ _: 

The stress-energy tensor component T'^^ is given by 



1 f Ar 



TT \ ab 



1 e' 



Substituting (|) and (|) in (0) we obtain 

A = 



a b 

Rearranging terms in equation ( p_0[ ) and using (|9|) we obtain 

1 



2r" 



2a'r' 



+ Ar^r' 



Integrating and simplifying the previous equation yields 

2 

-I -1- 
a 



2M 



(7) 



(8) 



(9) 



(10) 



(11) 



(12) 



where 2M is an integration constant and a 

2 



-A/3. The integration constant 2M can be evaluated using the black 



hole topology condition 



0, 



2M = r+-\ h 



(13) 



This is the known relation between the ADM mass of the Rcissncr-Nordstrom-anti-de Sitter black hole and its event 
horizon radius. 

Substituting (|l^ in (|) yields 



^ 2 



2ir 



dr 



dy -2- 



- 2- 



- 2- 



dr 



dy 



brr' 



-^subtr 
1 



^2 

' A' ■ 

ab 



dr 



b'r^ 



'subtr ; 



(14) 



y=o 
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where we have used the topology conditions given in section || and to evaluate the term in Ar , we used (|), since this 
term is identical to Im- 

The first term after the second equality in equation can be evaluated by integrating and substituting equations 
and ([ij). The respective third term is integrated and using the topology conditions gives — 7rr^. Following this 
procedure, we obtain 

I* = -f3rB f {rs; r+, e, a) - ^ /? e - 7rr^ - /subtr , (15) 

where the inverse temperature at the boundary (3 is given by the proper length of the time coordinate at the boundary, 
/3 = T-i = /p''6(l)dr = 2 7rfo(l) and 



f{rB;r+,e,a) = \ll-^-^^-a^^ + ^ + a^rl . (16) 

rB r+TB TB 



Adding (|6|) and (|15|), yields the reduced action 

r ^ -l3rB f [tb; r+, e, a) ~ [iecj)- TTr\ - I^nhtr ■ (17) 

The term /subtr is of the form /J-Esubtr, where -Bsubtr is a constant that does not depend on [3 or 0, since /subtr is an 
arbitrary term that can be used to fix the zero of the energy but cannot affect other thermodynamical variables 
For convenience, we use for the zero of the energy 

Sads = E{r+ = 0, e = 0) = , (18) 

where -Bads is the thermal energy of anti-de Sitter spacetime. 

To evaluate /subtr, we compute the thermal energy of the Reissner-Nordstrom-anti-de Sitter black hole from (|l 
and use condition dlq). The thermal energy is given by pq| 



E = F + f3 



/dF\ _f^\ 



= -rB f [rB; r+, e, a) - E'subtr (19) 

where F is the grand canonical potential and we have used (^. Although the reduced action /* is not the classical 
action (therefore we cannot write /* = (3F), the energy has the form given in (pj|). This is because the classical action 

I corresponds to the minimum of the reduced action and therefore the equalities ( §i ) = ( ) and 

(if) ~ (^^) hold. However, r+ and e in (19) are not free parameters, they depend on the boundary 

conditions (i.e. on the values of /?, (j) and r^) and on the cosmological constant. The functions r_|_ = r+(/3, 0, r^, a) 

and e — e{f3, 4>, tb, a) are obtained from the equilibrium conditions = and = as will be seen later. 
Using equation (18) on (|l^), yields 

-Esubtr = -rB fo [rB] a) (20) 



where fo (rB] a) = f {tb; 0, 0, a) = ^/T+oP^. 

Substituting (^^ in (p7|), we finally obtain the reduced action for the Rcissncr-Nordstrom-anti-de Sitter black hole 

r = [3rB [fo {rB; a) - f {tb; r+, e, a)] - /3e^- 7rr+ . (21) 

Similarly substituting (^0|) in (|l9|), we obtain its thermal energy 

E = rB [/o (rs; a) ~ f (r^; r+, e, a)] . (22) 

The mean value of the charge is 

fdF\ 1 ( di\ 

Q = -hjT ^-^Ut =e. (23) 
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The entropy is obtained from 



where equation was used. Since A^/A — vrr^, where is the area of the event horizon, this is the usual 
Hawking-Bekenstein entropy [ p^|j2^ . 

As mentioned above, the event horizon radius r+ and electric charge e of the black hole for the given boundary 
conditions, i.e. /3, (j) and vb, are obtained by evaluating the locally stable stationary points of the reduced action 
with respect to r+ and e Effectively, once the values of (3, 4> and tb are held fixed by the boundary conditions, 
then the reduced action is a function of only tj^ and e, i.e. /* = I*{rj^, e). The local stability conditions are then (i) 
V/* = and (m) the Hessian matrix is positive definite. The latter condition corresponds to a condition of dynamical 
as well as thermodynamical stability and will be discussed in section 0. We will start by investigating the first 
condition. 

The condition of stationarity V/* = gives 

(3 = — firg; r+, e, a) , (25) 

K 

2 2 I o 2 4 

where k = — — ^-^ is the surface gravity of the horizon, and 



f{rB;r+,e,a) . (26) 

r+ tbJ 

These are the inverse Hawking temperature and the difference in the electrostatic potential between r+ and tb 
blueshifted from infinity to rs, respectively. 

Inverting these two equations, r+ and e are obtained as functions of the boundary conditions and the cosmological 
constant. We define the more convenient variables 

x = — , q= — , /3 = , a =a Tb ■ (27) 

tb tb 47r rs 



In these variables (|25|) and (26) are written 



/3 = 2;VT~^\/l-^+s2(l + x + a;2) (^l_l_ + 3a2^2 j ^ ^^^^^ 



(t)=-Vl^[l- — +o?{l + x + x'^)\ . (29) 

X \ X J 

To invert this equations, we start by inverting equation ( p9|) to obtain q, 

^x^^^ [l + a^{l + x + x^)] {1-x + x 0^)-! . (30) 



-1/2 



Substituting ( pOD in (|2^) and taking its square, we obtain a 7th degree equation in x with a double root x = 1 which 
is not a solution of the initial equations. Getting rid of that root, we obtain a 5th degree equation 

+ |^9S2;^2 _ ^2 _ -^g-2^2^2 _^g_2■^2^4^| ^4 ^ (l - 0^) a;5 ^ . (31) 

However not every solution of this equation corresponds to a physical solution of a black hole. This is because the 
radius of the event horizon of the Reissner-Nordstrom-anti-de Sitter must obey the following condition 

rl-e^ + 3aV4 > o . (32) 



' + 



Where the equality defines the extremal Reissner-Nordstrom-anti-de Sitter black hole. 



5 



Comparing (32) with (^5|), yields that (3 is real and positive. Comparing it with ( p6[ ) we obtain the following 
condition 



1 + SaV^ 
l + a2(l + 2r+ rB + irl) 



< . , ..n ... 1 . o.2^ • (33) 



In the coordinates given in (|27|), the inequality becomes 



1 + m^x^ 

l + a2(l + 2a; + 3x2) 



This is the condition that the solutions of equation (^1|) must obey in order to represent physical black hole solutions. 

Equation (|3^) has no known analytical solutions. However its solutions can be numerically computed and presented 
in graphics. This will be done in the next section. 



IV. ANALYSIS OF THE BLACK HOLE SOLUTIONS 

In this section we present in graphics an analysis of the solutions of equation ( |3T| ) that obey condition (p^. This 
analysis is done in two steps: (i) first, we analyze figures |l] to ^, that present the solutions x as functions of P and 0, 
for values a = 0, 0.5, 1 and 5; (ii) afterwards we show in figures ||to I the regions with zero, one and two solutions 
in the space spanned hy (f> xa for fixed values of /3. 

(i) Analysis of figures |^ to ^: 

a — 0: The solutions for a = are presented in figure 111. These are obviously identical to the solutions of the 
Reissner- Nordstrom black hole, see (|ll|. We can see in figure^ that for fixed (p there is a maximum of /3, l3^^^{(f>), so 
that for P > /3,iiax(^^) there are no solutions. For /3 < /3,„ax(<^) can have two or only one solution depending on 
the precise value of (3. For = (Schwarzschild) one has always two solutions for (3 < /3max(0) (/^max(O) = 2/\/27). 
In the limiting case /? — > (i.e. rsT —> oo), one finds there is a solution with x = r-^^/rB 1 (as will be seen in 
section^ this is the stable solution), see js). For = 0, still, and /3 > /?niax(0) there are no solutions (see comments 
at the end of this section). For < < l/-\/3 one has one or two solutions up to (3^^^{(j)), whereas for (3 > Pj^axi't') 
there are as well no solutions. Finally for (j) > l/\/3 there is only one solution (again for (3 < I3^^^{(j))) corresponding 
to the unstable branch as will be seen in section ^ Note that, for a = 0, condition ( ^4|) implies that the electrostatic 
potential has a maximum at (/>max — 1- Notice also that in the limit /3 — > oo (T — > 0), the curves in figure |l| tend to 
the curve 0=1, which corresponds to the extremal Reissner-Nordstrom black hole (r-|- = e). 

The cases a ^ 0, which are now going to be analyzed require the following prior analysis. 

As in the case a = 0, for a 7^ there are solutions at T = (/3 = 00), that correspond to the extremal black holes. 



This can be analytically verified by replacing j3 — 00 isi equation (31), from where we obtain the equation 



1 - 0^ - a^0^ - 2S^0^x + 3a^(l - (j)^)x^ = . (35) 



Notice this is the equation one obtains taking the equality in condition 
of cxtrcmality of the Reissncr-Nordstrom-anti-de Sitter black hole, whie 
that only the extremal black holes have zero temperature. 

Equation (BS) has at least one solution that verifies < x < 1 if 



In fact it corresponds to the condition 
1 is in agreement with the well known fact 



2 1 1+3S2 ,n 1 

a < - and ^ < ch'^ < 



6 i + es^ 

^ <-o? <^ and ' ' < 0^ < 



1 _2 2 , 3^2+6 - VScr' + l2o? 



^^>l and ^ + ^^^!<02<i±i^. (36) 



For these values of and a, the curves for fixed presented in the figures reach infinite (3. Furthermore, equation 
(p5|) has two solutions if 
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1/6 < < 2/3 and ^ < cf) < ^j? , 

Aa +6 1 + 6a 



> 2/3 and ^ < c^^ < ^ . (37) 

a — 0.5: The solutions for a — 0.5 are presented in figure |^. This figure presents the same properties mentioned 
previously for the case a ^ Q. Comparing |l| and H, we verify that the maximum value of /3, /?niax('/')i for which there 
are solutions is increasing, i.e. there are solutions for slightly lower values of the temperature at a = 0.5 than at 



a = 0. Furthermore there are solution for infinite (3 in the interval 0.83 < < 0.89, see (36), this can be seen using 
the curve /3 = 9 in figure ^, since it corresponds to a good approximation of infinite 
Condition (|3|) implies, for < 2/3, 



This condition corresponds to the upper-limit of the interval given in (|36D, which in this case is </) ~ 0.89. Notice that 

if in equation ( |3l] ) we do a; = 0, we obtain precisely (j) — 1 / \/ 1 + a ^ , as can be seen in figure |[ 

a — 1: Figure |^ plots the solutions for a = 1. We can see that this figure presents similar properties to the ones 
obtained for lower values of a. Using (^6|), we can see that for cj) < 0.66 there are solutions only for /? < /?niax(</'); 
where /3max('/') finite and depends on (j). On the contrary, the curves for higher values of (j) reach infinite (3. In 
particular, for 0.66 < </> < 0.71 (see equation (^) there are two solutions at low temperatures (i.e., high /3) as can 
be seen in figure ||, since the curve /? = 9 is representative of the curves with high (3. It can be seen that for /3 = 9 
there are two solutions for (f) < (j)o, where (po = , ^ _ ~ 0.71 is the value of where a; = for every /3. There is one 

solution for 0.71 < (/i < 0.76, see figure ||, where the upper-fimit is imposed by condition (^. In fact, condition ( |3^ ) 
imposes, for < 2/3, 



^ . (39) 



Notice this is the upper- limit of the interval given in (|3q). 

a = 5: Figure |^ presents the solutions for a = 5. Using (|3^), we can see that for ^ < 0.19, there are solutions 
only for /3 < /3i„ax(0)- For 0.19 < ^ 0.71 there is one solution for high values of /3 (consider the curve /? = 9). For 
infinite /3 this region is 0.195 < < 0.196, see (|37|). In figure ^ we can also see that for 0.19 < </> < 0.71 there is one 



solution, where the upper-limit is given by (39) 



For higher values of a there are not new types of solutions and therefore it is not necessary to pursue our analysis. 
ill) Analysis of figures ^ to ||: 

In order to clarify the disposition of the number of solutions for given values of a, ^, and /J, we present, in the space 
spanned by a and for fixed /?, the regions with zero, one and two solutions. We do this for eleven different values 
of /3 = 0, 0.3, 2/\/27 ~ 0.38, 1, 00, see figures || to ^ respectively. In this figures one can see the evolution of the 
number of solutions as /3 increases. To present all possible values of ae[0, 00 [, we use in these figures the parameter a 
instead of a, a is defined by 

2 _ 

a — — arctan a , (40) 

TT 

such that < a < 1. It is this variable that appears in the ordinate axis in figures H to ^. 

Due to condition (^4|) there are no physically possible solutions on the right-hand side of figures ^ to |[ 

An important value of /3, first studied by York [|| in connection with the Schwarzschild black hole (0 = a = 0), is 

/3 = 2/V27 ~ 0.38, i.e. /3 = tb- For higher values of lower values of the temperature, there are no black hole 

solutions. This is a quantum effect and following York |^ can be understood as follows. One can associate a Compton 
type wavelength A to the energy ksT of the thermal particles, by A = or in Planck units A = 1/T = p. If this 
Compton wavelength is much larger than the radius rg of the cavity (or more specifically A > rg ~ 4.8 r^) then 
the thermal particles cannot be confined within the cavity and do not collapse to form a black hole. 

By analyzing figures |^ and ^ we can see that for nonzero cosmological constant (a 7^ 0) this phenomenon starts at 
even lower P (higher T). Indeed using equation ( ^8|) (with (j) — 0, i.e. q — 0) we can show that to first order in 
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{a^ rg << 1), York's criterion for no black hole solutions becomes 

From equation (^) we infer that the role of the negative cosmological constant (A = — 3a^) is to produce an effective 
cavity radius r^g = (l — jgCt^ r^) smaller than rg. Thus for a given temperature, it is more difficult to confine 
the thermal particles, and harder to form black holes, in accord with the idea that a negative cosmological constant 
shrinks space. 

If we extend our previous first order analysis to include the electrostatic potential (jj, we obtain 

We can see that the electrostatic potential has the opposite effect of the cosmological constant (see for example figures 
II and |. 



V. STABILITY 



A. Local stability 



As mentioned before there is a second condition of local stability that has not yet been investigated. We will follow 
the same procedure as given in p^ . This is the condition that the Hessian matrix of the reduced action be positive 
definite. For convenience in this analysis we will use the variable S — 7rr^ instead The Hessian matrix of /*(S', e) 
is 



eS 



The matrix is positive definite if its pivots are all positive. The pivots of ( ^3| ) are 

/* 

,ee 

and 

det(/*,) 



(43) 



(44) 



(45) 



The first condition of local stability, V/* = yields 









de ) 


^ de 






^dE* 




V dS 



/30 = 



1 = 



dE 



s 

dE 
95 



These are well known thermodynamical relations [ [TgI . 

From (^), we compute the second derivatives of /* in the stationary points of /* 



Q2p 




d^E* 




cq 


= ^ de^ 






d^E* 


de dS 


oq 


~ ^ dedS 






d^E* 




eq 


~ ^ ds^ 



cq 



^ ' de 



= /3 



90 
'dS 

1 (di 

[3 \dS 



1 



dp_ 

de 



(46) 



(47) 



where eq means quantities evaluated at equilibrium, i.e. at the stationary points of the reduced action /*. 
The first pivot (H) is simply the first of these equations. The second pivot (|4|) is easily obtained from (|47 
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detl* 



dl3_ 
dS 



dl3_ 
dS 



dp 



f3 \dS 



1 /dp 

13 \de 

d(3_ 
de 

1 

C(f>,rB 



dS 



de 



de_ 
dS 



(48) 



where C^^rs is the heat capacity at constant <j) and tb- Notice that in aU this calculation we have implicitly held tb 
constant, since in the grand canonical ensemble the dimension of the system is held constant. 
Imposing positive pivots for local stability yields 



(49) 



These conditions are identical to the classical thermodynamical stability conditions jl^ . Therefore one can conclude 
that the dynamical stability conditions given in ( ^ ) and ( ^9|) are identical to the thermodynamical stability conditions 

For the Reissner-Nordstrom-anti-de Sitter black hole, the pivots obtained from (E^ and (Bs) are 



9^ 
de 



S,rB 



1 






_'[± _ 


e^ 


r+ 


TB J 


rB 


r+TB 


e2 


it 






!± _ 


rB 




TB 


rB 



TB 



Tb 
rB 



2 2 

a rB 



1/2 



2 2 

a rB 



3/2 



(50) 



which is positive, and 



^,rB = 47rr+^(r - r+){r+^ - + 3aV+'*)[l + a^{rB^ + r+rs + r+^)] x 
X {e^ + r+3[(6QV+2 - 2){rB + a^rB^) + 3r+ + 2a^r+^ + 3aV+5] + 



+ 2eV+[-2r+ + rB{l + a^{rB^ ~ 2r+ri 



2r. 



(51) 



The numerator of C^^rs positive, therefore the condition C^ ^b > is verified if the denominator in ( [5l| ) is positive. 
Using (p7|) and (pQ), we obtain the following condition of stability for the solutions of equation (|3l|) 



[-2 (1 + 3.T + 65^(1 + 0^) + 2s2x^ + 3a^x^] (1 - a; + x^^)^ - 

- [2 02(1 + a2(l + X + x^)){-l + 2x + a^{-l + 2x + 2x'^))\ (1 - a; + x<f) + 
+ (j)'^x{\ + a^{l + x + x^)f > 0. 



(52) 



By numerical computation, we can verify the solutions that obey this condition. We have found that the lower branches 
of the curves presented in figures |l| to ^ correspond to unstable solutions, while the upper branches correspond to 
stable solutions. Therefore we can say that in general only the solutions with the higher value of x, that is with higher 
event horizon radius, are stable. 

In more detail we can distinguish three cases: (i) for low a (see figures Q and and for values of (3 and (j) for which 
there are one solution, it corresponds to a lower branch and therefore this solutions is unstable; (ii) for high values 
of a (see figure^), whenever there is only one solution, it corresponds to an upper branch and therefore this solution 
is stable; (iii) whenever there are two solutions, for given values of (3,(f>,a, the smaller one is an instanton (i.e. it is 
an unstable solution and dominates the semi-classical evaluation of the rate of nucleation of black holes ||2^) and the 
one with larger event horizon radius is a stable black hole. 

These three cases can easily be distinguished in figures || to ^. In these figures there are in general two separated 
regions with one solution. The region with lower values of a, i.e. lower values of a, corresponds to case (z) and these 
are unstable solutions. The region with higher values of a corresponds to case (ii) and these are stable solutions. 
Obviously, the region with two solutions in each figure || to ||, corresponds to case {in), i.e. one of those solution is 
stable, the one with higher value of x, and the other unstable. 
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B. Global stability 



The stable solutions computed in the previous subsection are not necessarily global minimum of the action. In this 
case they will not dominate the partition function and the zero-loop approximation cannot be considered accurate 

10- . . . . . 

As the reduced action given in equation (Ell) grows without bound in the directions where tj^ or e tend to infinity, 

the global minimum must be found either at the local minimum or at r_|_ = e = 0. The action tends to zero as r_|_ and 
e tend to zero, therefore the condition of global stability of the stable solutions is that the classical action / must be 
negative. 

If the classical action is positive, the partition function is dominated by points near the origin. But these points 
do not correspond to a black hole in thermal equilibrium. In this case, the black hole, that corresponds to the stable 
solution of the reduced action, is metastable. 

We verified numerically which boundary conditions (given by the values of /?, 4> and a) correspond to globally stable 
black holes, i.e. to solutions that dominate the partition function. We can also see, using a simple argument, that 
for /3 < ^ all locally stable solutions are also globally stable. Indeed, York [|j has shown that for the Schwarzschild 
black hole the condition of global stability \s (3 < a-nd since the classical action / is a decreasing function of (f) and 
a, this condition is still a bound for globally stability for all cf) and a. However for (5 > there is always a certain 
region oi (p x a, for which the locally stable solutions do not correspond to global minima of the action. We show 
these regions in figures ^ to where the regions for which the solutions do not correspond to a global minimum of the 
action, i.e. the metastable solutions, are shaded. In particular figures ^ and |9| do not present a shaded region because 
all the stable solutions are dominant in the limit /3 ^ as said above and for (3 ^ oo the region with metastable 
solutions is too thin to be presented in graphic. 



VI. THE Rb ^oo LIMIT AND THE HAWKING-PAGE SOLUTIONS 

One can study the case where the boundary goes to infinity. There are two different ways for taking this limit: («) 
fixing the horizon radius and the charge e of the black hole; (m) fixing the boundary conditions, i.e. fixing p and 
(j). In this section we will study only the pure Reissner-Nordstrom-anti-de Sitter cases, i.e. the cases with a,(j) ^ 0. 
Other cases, with either a or (f) equal to zero are considered in the next section. 

We will start by studying the first case: (i) Fixing the black hole solution, i.e. fixing r+ and e and taking the 
limit vb oo, we obtain from equations (E3) and (pq) that the temperature T — (3^^ and electrostatic potential (f> 



go to zero as T inar^ ^^"^ ~ ^B^' respectively. In this case the thermal energy goes to zero as 

^b) ' where Af is the mass of the black hole given in (|l3|). As also previously found in |l2| , the 
thermal energy at infinity is not equal to the ADM mass of the black hole, due to the the fact that the spacetime is 
not asymptotically flat. Note that in (l^Jl^ it is found that E = M due to a different definition of the temperature 
of the ensemble. 

To determine the stability of the black hole solutions in this limit, we compute the heat capacity, from equation 
( pl| ) and obtain 

^<t> — 5 2 ~o 2~T ■ ^^'"^"i 

The numerator in ( p3[ ) is necessarily positive due to condition (p2[), that r-|_ must obey to represent the event horizon 
radius. Therefore, the stability condition for these solutions is simply — + Sa^rl > 0. Notice that equation 
( ^ ) is a generalization to e 7^ of the heat capacity found by Hawking and Page in |p^ . Computing the action in 
this limit yields 



{r\ — — 



3 rl 



(54) 



This is precisely Hawking-Page action [^7|jT^] , which here we have recovered from York's formalism taking the appro- 
priate limit. 



From equation (54), we can verify the global stability of the black hole by imposing that / < 0. Therefore we 
conclude that the solution given by (r+ , e) corresponds to a globally stable black hole if conditions — — < 

and ([33) are both verified. These conditions are similar to those found in iflSl . 
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We conclude that taking different boundary conditions, i.e. cfioosing to fix the boundary conditions in the boundary 
at infinity (as done here) or in the region where spacetime is flat like Hawking and Page ||l7[] , yields a different value 
for the energy only. Furthermore this difference is so that all the other physical quantities (like the action, entropy, 
mean- value of the charge and heat capacity) remain the same. 

On the other hand, we can take a different limit: (ii) we can fix the boundary conditions, i.e. (3 and 4>, when taking 
the limit. This is done by recovering the variables 13 and a in equation (|3^), using (p7|). Taking the limit tb ^ 00, 

we obtain the equation (0^ {l + 2x + ix^) — 3x^) — a;^ (1 + x + x^) {1 — x + (fP' x) — 0. This is a 5th degree 

equation in x. For fixed /3, <p and a and taking the limit rs —>■ 00, we obtain a solution r+ that tends to infinity like 
r+ ~ c+r^, where c+ is a constant that depends only on (3, (f) and a. Considering now equation (^0|), we can see that 
the charge e goes to infinity as e '--^ CeTg, where again Cg is a constant that depends only on (3, (f) and a. Therefore 
the entropy ( P4[ ) and the mean- value of the charge ( |23|) b oth go to infinity as r^. In this limit the action, the thermal 
energy and the heat capacity, given in ( pl| ) , (^2|) and (|5l|) respectively, also diverge as . The heat capacity is always 
positive, which means these solutions are stable. The action is always positive, therefore the solutions are not globally 
stable and represent metastable black holes. 

VII. COMMENTS ON SPECIAL CASES 

Several black holes may be considered as special cases of the Reissner-Nordstrom-anti-de Sitter black hole. 

(I) putting (p — and A = 0, we obtain the Schwarzschild black hole studied in ||]. There are two solutions for 

f3 < (3^^^ — This solutions can be computed analytically, since equation ( pT| ) becomes a 3rd degree 

equation for A = 0. Only the solution with higher event horizon radius, i.e. higher mass, is stable. 

In the limit rs — > 00 the unstable solution, i.e. the solution with lower value of the horizon radius, goes to 
r+ — l/(47rr+), while the stable solution goes to infinity as r+ ^ H]. 

(II) putting A = we obtain the Reissner-Nordstrom black hole. This has been studied in There are one or two 

solution for < f3^^^. These as stated above can be computed analytically. Once again only the solution with 
higher event horizon radius, when it exists, is stable. 

In the limit — > 00, the black hole horizon radius is r_|_ = ^ (l — </)^) and the charge is given by e = ^t^. 
The thermal energy of these solutions is equal to their ADM mass E = M. The heat capacity is negative 
— —2-Kr\. therefore the solutions are unstable. 

(III) putting (/) = we obtain the Schwarzschild-anti-de Sitter. This black hole has been studied before in |]l7| , |l2| . 
This black hole has two solutions for (3 < /Jmax, and again only the one with higher event horizon radius is 
stable. 

The limit — s- cx) can be taken in two ways: (i) fixing the temperature and the cosmological constant, there 
are one unstable solution that tends to zero as r+ ~ "S^^b^ ^^'^ o^*^ stable solution that tends to infinity as 

r+ ~ crB, where c is the solution of equation + ^^^^ ~ I — 0; (ii) fixing the horizon radius r-|_, the 
temperature goes to zero as T ~ ^4^ar+^ ^b^^ these solutions are stable if 3a^r^ > 1, sec Jl2| . 

(IV) the extremal cases require special care p3[-p5[l and were not studied in any detail in this paper. 

VIII. CONCLUSIONS 

We have studied the thermodynamics of the Reissner-Nordstrom-anti-de Sitter black hole in York's formalism. In 
the grand canonical ensemble where the temperature and the electrostatic potential are fixed at a boundary with 
finite radius, we have found one or two black hole solutions depending on the boundary conditions and the value of 
the cosmological constant. In general when there are two solutions, one is stable, the one with larger event horizon 
radius, and the other is an instanton. On the other hand, the cases with a single solution can correspond either to a 
stable or unstable solution. We have found that both high values of the cosmological constant and low temperatures 
favor stable solutions. 
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FIG. 1. Solutions of equation (^l|) for a = (Reissner-Nordstrom) as a function of the electrostatic potential at the boundary 
(j> for fixed values of /3 = 0.1, 0.3, 0.6, 0.9, 3, 9. The stable solutions correspond to the upper branch of the curves. This means 
that when there are 2 solutions for given values of P and 0, only the solution with higher value of x is stable. 



/3=0.1 

^=0.3 
---- 0=0.6 




0=0.9 
^=3 
0=9 



0.5 1 

FIG. 2. Solutions of equation (^l|) for a = 0.5 as a function of the electrostatic potential at the boundary (f> for fixed values 
oi (3 — 0.1,0.3,0.6,0.9,3,9. Notice that <j!> < .89, as imposed by condition (p^. The stable solutions correspond to the upper 
branch of the curves. This means that when there are 2 solutions for given values of /3 and (p, only the solution with higher 
value of X is stable. 
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FIG. 3. Solutions of equation ( ^ ) for q = 1 as a function of the electrostatic potential at the boundary (p for fixed values 
of /3 = 0.1, 0.3, 0.6, 0.9, 3, 9. The maximum value of </> for which there are solutions, t/imax — 0.76 is imposed by condition (^^. 
The stable solutions correspond to the upper branch of the curves. This means that when there are 2 solutions for given values 
of /3 and (j), only the solution with higher value of x is stable. 
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/3=0.1 

^=0.3 

____ P=0.6 
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FIG. 4. Solutions of equation ( ^ ) for a = 5 as a function of the electrostatic potential at the boundary (j) for fixed values 
oi P — 0.1, 0.3, 0.6, 0.9, 3, 9. The maximum value of (p for which there are solutions, (^max — 0.71, is imposed by condition (^). 
The stable solutions correspond to the upper branch of the curves. This means that when there are 2 solutions for given values 
of P and (/), only the solution with higher value of x is stable. 




FIG. 5. Number of solutions of equation ( |31| ) with /9 — ^ 0, in the space x a, where a is given by equation (po| ). There is 
one black hole solution in the confined region and also for a — 1 (i.e. infinite cosmological constant), for (f) < VO.S. There are 
two solutions for 6 = 0, i.e. for the Schwarzschild-anti-de Sitter black hole. 




FIG. 6. Number of solutions of equation (^l|) with /3 — 0.3, in the space x a, where a is given by equation ( ^ ) - means 
that there are zero solutions in this region, i.e. there are no solutions of black holes in thermodynamical equilibrium for this 
set of values of /?, cj) and a. 1 - means that there is one solution. 2 - means that there are two solutions. In th e shaded region 
the stable solutions are not globally stable and therefore represent metastable black holes (see subsection VB). 
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FIG. 7. Number of solutions of equation with /3 = ~ 0.38, in the space <f> x a, where a is given by equation ( |4c| ) 
(see caption of figure m). 




FIG. 8. Number of solutions of equation with /3 = 1, in the space (j) x a, where a is given by equation ( pc| ) (see caption 
of figure 0) . 




FIG. 9. Number of solutions of equation (^ with /3 = cxa, in the space x a, where a is given by equation (^). Notice that 
for /3 — + cx), equation (js^) becomes ( ^ ) which corresponds to the extremal Reissner-Nordstrom-anti-de Sitter black hole, see 
discussion following equation (ISq) (see caption of figure H). 
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